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We combine a conventional harmonic analysis of vibrations in a one-atomic 
model glass of soft spheres with a Voronoi-Delaunay geometrical analysis of 
the structure. "Structure potentials" (tetragonality, sphericity or perfectness) 
are introduced to describe the shape of the local atomic configurations (Delau- 
nay simplices) as function of the atomic coordinates. Apart from the highest 
and lowest frequencies the amplitude weighted "structure potential" varies 
only little with frequency. The movement of atoms in soft modes causes tran- 
sitions between different "perfect" realizations of local structure. As for the 
potential energy a dynamic matrix can be defined for the "structure poten- 
tial". Its expectation value with respect to the vibrational modes increases 
nearly linearly with frequency and shows a clear indication of the boson peak. 
iq \ The structure eigenvectors of this dynamical matrix are strongly correlated 

to the vibrational ones. Four subgroups of modes can be distinguished. 
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The thermodynamic properties of glasses at low temperatures differ from those of 
the corresponding crystals^. At low temperatures the specific heat is strongly enhanced 
compared to the Debye contribution stemming from the sound waves. The excitations 
underlying this enhancement have been shown to be two-level systems below T « 1 K and 
nearly harmonic vibrations above. The vibrational density of state, Z{y), plotted as Z[y)jv 2 
rS \ has a maximum, typically near 1 THz, the boson peak. 

This low temperature / low frequency behavior can be described by the soft potential 
model.0'0 In this model one assumes that one common type of structural unit is responsible 
for the excess excitations. One introduces an effective potential describing the motion of 
this unit. Depending on the parameters this potential is a single well or a double well. In 
the first case it describes a low frequency localized vibration and in the second tunneling 
through the barrier (two level systems) or relaxation over the barrier. For low energies 
one can give a general form for the distribution of the parameters describing the effective 
potentials. Fitting this model to the experimental data, one finds that 20-100 atoms or 
molecular units move collectively in the tunneling and in localized vibrations.Q'Q It should be 



emphasized that the concept of low frequency localized vibrations is an idealization. These 
modes will always interact with the sound waves of similar frequency and, therefore, also 
among each other. This delocalises the modes and they are only quasi- local or resonant. 
Due to level repulsion, for sufficiently high densities of these modes, the interaction will 
change their density of states from Z{y) oc z/ 4 to Z(y) oc v thus creating the boson peakB 
Such a model does, however, not say anything about the physical nature of the localized 
modes or their origin in different types of glass. 

The problem of local dynamics in the amorphous state is closely connected with the 
problem of the so-called medium-range order in glassesQ'B. Recently it has been showntiI2l 
that a computer model of amorphous argon has a heterogeneous structure containing regions 
of more "perfect" or "imperfect" atomic arrangements on a nanometer scale. In the regions of 
perfect structure the elementary packings of four neighboring atoms (the Delaunay simplices) 
are close to either regular tetrahedra or quart-octahedraa, i.e. quarters of regular octahedra. 
In the regions of imperfect structure the local configurations of the neighboring atoms differ 
markedly from these ideal shapes. A partial spectrum of the vibrational states of the atoms 
in the regions of more "imperfect" structure displays an excess for low-frequency modes." 

Quasi-localized low-frequency vibrations have been observed in computer simulations of 
the soft sphere glass (SSG)EJ and of numerous other materials, such as e.g. Si02^, Se» 
in Ni-Zio and Pd-SicJ, in amorphous iceE^I and in amorphous and quasi-crystalline Al-Zn- 
MgO. It was shown that these modes are centered at atoms whose structural surrounding 
differs substantially from the average. It has been established that the directions of the 
eigenvectors of soft vibrations strongly correlate with those of the relaxation jumps at low 
temperatures.^^ 

One hypothesis on the origin of the soft mode is that the most active atoms oscillate be- 
tween neighboring minima of the potential energy formed by a cage of surrounding atoms£2l 
These minima correspond to some more "perfect" local arrangements of the atoms. The 
coupling to the rest of the material changes this double well system to a soft single well 
one. One example for such a situation is the interstitial atom in an fee metal. A medium 
sized interstitial occupies the octahedral site. Increasing the size of the interstitial atom the 
octahedral site becomes unstable and the interstitial moves to an off-center position. The 
impeding instability is indicated by low lying resonance vibrationsciEHl. The instability in 
this example is caused by a local compression which causes the simultaneous occurrence 
of high frequency localized vibrations. In the glass the modes are more extended typically 
string like groups of some twenty atomsc3. Instead of the single interstitial atom one has to 
take a group of atoms and due to the laking symmetry the energy minima will be shifted 
relative to each other. Keeping this in mind the underlying mechanism can still be true. 
The simultaneous occurrence of low and high frequency localized modes centered on one 
atom has indeed been observedllil. 

In the present paper we want to verify and concretize this notion for the SSG. For this 
purpose we combine the harmonic analysis of Ref. |ll] with the Voronoi-Delaunay geomet- 
rical description of the local structure used in Ref. |TU[ First we shift the atoms of the 
model along the eigenvector of a low frequency quasi-localized normal mode and observe 
the changes in the local atomic arrangements, caused by the shifting. This allows us to 
visualize the specific transformations of local structure which accompany the movement of 
atoms in the soft vibrations. In a next step we calculate the atomic perfectness weighted 



with the squared amplitudes of the vibrational modes. This quantity varies only weakly 
with frequency. Considering that the vibrations are connected with changes in the geome- 
try, we introduce a "structural dynamical matrix" . We will show in the following that there 
is a strong correlation between the "structural eigenvectors" and their vibrational counter- 
parts. This correlation divides the vibrations, as regards structure changes, into separate 
classes: longitudinal and transverse extended, high frequency localized and low frequency 
quasi-localized modes. 

II. THE SOFT SPHERE GLASS 

We use 55 glassy configurations of 500 atoms each, interacting via a soft sphere pair 
potential 

u{r) = e(^)\A^-f + B. (1) 

To simplify the simulation the potential is cut off at r/a = 3.0 and shifted by a polynomial 
with A = 2.54 x 10 _5 e and B = —3.43 x 10 _3 e. The calculations are done with a fixed atomic 
density, p/a 3 = 1 and periodic boundary conditions. The configurations were obtained by a 
quench from the liquid to T = K. From the pair correlation one finds a nearest neighbor 
distance of around 1.1c For more details see Ref. [TTJ. 

The inverse sixth-power potential is a well-studied theoretical model that mimics many 
of the structural and thermodynamic properties of bcc forming melts including the existence, 
in its bcc crystal form, of very soft shear modesEa In the glassy structure one finds a boson 
peak with a maximum near v = 0.1 (e/ mer 2 ) 1 / 2 extending to about v = 0.4(e/mcr 2 ) 1 / 2 . The 
enhancement of the vibrational density of states over the Debye value is by a factor of 2.5.E± 

As before the frequencies and eigenvectors of normal vibrations are calculated by the di- 
agonalisation of the force constant matrix. Imaginary frequencies are absent in the spectrum 
because the system is in an absolute local minimum of potential energy. For the given number 
of atoms the minimal q- value for sound waves is q m [ n = 0.79a -1 giving minimal frequencies 
of 0.18 and 0.62 (e/ma 2 ) 1 ^ 2 for the transverse and longitudinal sound waves, respectively. 
Resonant modes with frequencies well below 0.18 (e/ma 2 ) 1//2 will, therefore, be seen as low 
frequency localized modes. This is reflected in the participation ratios given in Ref. [TT| . 
One finds proper localized modes at frequencies v > 2(e/m<r 2 ) 1//2 and (quasi-) localized low 
frequency modes with v < 0.2(e/mo" 2 ) 1 ' 2 . The great majority of modes (0.2 < v < 2) 
extends over the system. These latter modes have been called diffusonsE3 due to their non- 
propagating character. Nevertheless for the SSG as for other systems it ispossible to extract 
via the dynamic structure factor some very broad "phonon dispersions" £3. 

The SSG was used in extensive studies of the influence of quench rate on the glass 
structureHH In these studies the Voronoi method was used to identify pentagonal rings 
which can be used as signature of icosahedral packing. 

III. VORONOI-DELAUNAY DESCRIPTION OF LOCAL STRUCTURE. 

By definition, the Voronoi polyhedron (VP) of an atom is that region of space which is 
closer to the given atom than to any other atom of the system. A dual system spanning 



space is formed by the Delaunay simplices (DS). These are tetrahedra formed by four atoms 
which lie on the surface of a sphere which does not contain any other atom. Both VP 
and DS fill the space of the system without gaps and overlaps. In our calculations we do 
a Voronoi-Delaunay tessellation of the glass configurations by the algorithm described in 
Ref. 0. 

It was found earlier that two main types of DS are predominant in mono-atomic 
glassesE5l'Ei namely DS similar to ideal tetrahedra and DS resembling a quarter of a regular 
octahedron (quart-octahedron). Following Ref. |32] we introduce as quantitative measure of 
tetragonality of a DS 
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where i and j designate the edges of the simplex, and I is the average edge-length. This 
measure was constructed to be zero for ideal tetrahedron and to increase with distortion. 
For computational reasons we slightly modify the previous measure of octagonality using: 
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In a perfect quart-octahedral DS one edge is \/2 times larger then the other edges. In 
Eq. dJ), the previously used measure^, it was originally assumed that the m-th edge is the 
longest. The modified expression (|3|) weights the six possible values O m in such a way that 
the smallest one dominates. The octagonality thus tends to zero when the DS is close to 
the perfect quart-octahedron. This weighting allows us to avoid the use of logical functions 
for the selection of the maximal edge and guarantees differentiability which is essential for 
our investigation. For the relevant low values of O, i.e simplices close to quart-octahedral 
shape, our expression reproduces the values of the original definition. 

The tetrahedral and quart-octahedral DS can be unified in one class of "perfect", or 
"ideal" simplicesa We measure the ideality of the DS shape by 
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S tends to zero when the simplex takes the shape of an ideal tetrahedron or quart- 
octahedron. Contrary to the expression proposed in Ref. ||] our measure is differentiable 
with respect to the atomic coordinates. The relative weights of tetrahedricity and octa- 



hedricity, T c = 0.016, O c = 0.033, are taken from Ref. [3^. The relation of the values T, O 
to the distortion of a DS can be also seen from the values To — 0.050 of the tetrahedricity 
of an ideal quart-octahedron and Or — 0.084 octahedricity of an ideal tetrahedron. 

Each atom in the glass is corner of approximately 24 DS. The structural environment of 
an individual atom can be characterized by the average ideality of over these DSu: 

J n DS 

'-'atom = /_, *->i (6) 

n DS U. 

where nos is the number of DS surrounding the atom. 

Another widely used measure of the atomic neighborhood is the sphericity of the Voronoi 
cell: 

1 F 3 

Sph = - 1. (7) 

y 36vrU 2 V ; 

Here F is the surface area of the VP, and V is its volume. This measure is minimal for a 
sphere , Sph = 0, and again increases with distortion. 

In analogy to the potential energy one can take the total tetrahedricity, ideality or 
sphericity to characterize the structure. We introduce an average "structure potential" by 

i n ds 
(T) = J7- E Ti (8) 

and analogously (S) and (Sph) where N^s is the number of DS in the system. Since 
dynamics is concerned with the motion of the atoms it is often more useful to average over 
the atomic quantities defined by Eq. || 

1 N 

\-L atomic/ i\t / . atomic- \J ) 

Both definition give similar values. 

In table [I] we compare the values of the three measures for an ideal fcc-structure an 
icosahedron and our glass. The values of the glassy structure clearly deviate from the ones 
of both the ideal configurations. It is, however, not possible to define unambiguously a 
nearness to either structure using these measures. 



IV. SOFT VIBRATIONS AND CHANGE OF THE LOCAL STRUCTURE 

Using the quantities defined above we will now illustrate for one example of a quasi- 
localized soft mode the relationship between softness and local geometry. In Fig. [l] (solid 
line) we show the average potential energy per atom as function of the displacement along 
a single soft eigenmode, i.e. one of the soft potentials which are described by the soft 



potential modeffl discussed in the introduction. The atoms are shifted along the direction 
of the 3A-dimensional eigenvector e as 

R n (x) = R™ + xe n . (10) 

Here Rq is the equilibrium position of atom n. For simplicity we have not normalized the 
amplitude x to an effective atomic amplitude as is usually done in the soft potential model. 
Fig. |l| corresponds to a very well localized soft mode with v = 0.04(e/mcr 2 ) 1//2 , effective mass 
13m and participation ratio 0.14. These values would guarantee a very narrow resonance in 
the infinite mediumBH. 

As mentioned in the introduction it has been speculated that the soft modes in glasses 
originate in some "soft" atomic configurations where, in the extreme case, the atoms are 
stabilized by the embedding matrix in a position lying between minima of the potential 
energy given by its near neighbors. In Fig. [I] we show by the dashed line the average 
potential energy of the 13 most active atoms of our soft mode, i.e. the atoms with the 
largest amplitude e n e n . This partial potential energy is indeed double- well shaped with 
minima at x m rs ±1.0 which corresponds to maximal displacements of individual atoms by 
|R ra — Rq| rj 0.27a from the equilibrium configuration. Note that at x = the selected 
atoms have somewhat smaller potential energy than the average. This reflects the reduced 
number of nearest neighbors reported earlier Jlil 

In order to understand the changes in the local structure as the atoms oscillate between 
the two partial minima of the potential energy, we look at the changes of the Delaunay 
tesselation caused by the displacements of the atoms in the normal mode. In Figs. |2] and 
|3] the five most active atoms of the soft mode are shown by the black spheres, and their 
geometrical neighbors by gray spheres. We consider two atoms as geometrical neighbors if 
they share a DS. The DS are visualized by line segments. 

In Fig. ^ we concentrate on DS with nearly ideal tetrahedral shape. For the sake of 
clarity only the most perfect tetrahedra (T < 0.003) are drawn. In equilibrium (x = 0) two 
tetrahedra are found which satisfy this condition (Fig. [|, b). After shifting in the "positive" 
direction (Fig. ||,a) one perfect tetrahedron has disappeared, but 6 new ones appeared. In 
particular, a 5-fold ring of perfect tetrahedra is created (seen sideways in the right half of 
the picture). This ring is known to be the densest possible packing of 7 equal spheres. After 
the displacement in the "negative" direction (Fig. [|,c) again one tetrahedron is lost, but 4 
new perfectly tetrahedral DS are gained. 

Similarly Fig. [^ demonstrates the appearance of new quart-octahedral DS in the neigh- 
borhood of the active atoms. Only DS with O < 0.008 are shown. For x = 1 one perfect 
octahedron, consisting of 3 active atoms and 3 of their neighbors, is observed. Together with 
the 5- fold ring of the perfect tetrahedra, it forms locally a pattern of perfect, non-crystalline 
structure which does not exist at x — 0. At x = — 1, a large dense cluster of 12 ideal 
quart-octahedra appears (Fig. [| c). It indicates several octahedral configurations in the 
neighborhood of the active atoms, although they cannot be seen clearly on the figure. 

In general, the number of perfect DS (tetrahedra or quart-octahedra) increases as the 
atoms are shifted from the equilibrium position to the local minima of the partial potential 
energy. This tendency is summarized in the double-well behavior of the average ideality, 
('S'atom), as function of the normal coordinate of the mode (Fig. f|, dashed line)t3. We 
remind that a lower values of S a t om means a more perfect atomic neighborhood. Note that 
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the minima of the partial ideality are situated approximately at the same values of x as 
the minima of the partial potential energy. In the equilibrium position (x = 0) the active 
atoms have relatively imperfect neighborhoods compared to the rest of the atoms. After the 
displacements this is considerably "improved" . 

The curve of the average perfectness (S at om(x)), averaged over all atoms (solid line) is 
almost flat at x = and resembles the behavior of the average potential per atom energy of 
the mode (Fig. [IJ, solid line). 

The double-well behavior of the partial potential energy and the ideality of the atomic 
environment is specific for a number of low-frequent vibrations. It becomes less pronounced 
as the frequency increases. However, we have not noticed a sharp transition between localized 
and delocalized low-frequency modes. Displacements of the atoms along the modes with 
medium and high frequencies also destroy DS of near perfect shape. 

The geometrical peculiarities of the random soft sphere packing play an important role 
for the localized low-frequent vibrations. These vibrations can be visualized as complex 
collective motions which organize the atoms in some "perfect" but non- crystalline arrange- 
ments. Although these arrangements consist of elements present in the crystalline structure 
(tetrahedra and quart-octahedra) these are connected to each other in a way which is incom- 
patible with rotational and translational symmetries of the crystal. Slight deviations of the 
shape of the DS make a variety of spatial arrangements possible which differ from the closed 
packed fee and hep ones. The pentagonal rings typical for locally icosahedric structure is 



one such example, compare e.g. Ref. 27. The lack of translational symmetry restricts these 



structures to ranges of a few interatomic distances. 



V. CORRELATION BETWEEN STRUCTURE AND VIBRATION 

We have seen that a low frequency quasi-localized vibration has a specific impact on 
the structure surrounding the most active atoms of the vibration. We will now investigate 
how far a general relationship between structural measures and dynamics can be seen. We 
will here concentrate on tetragonality, Eq. ^. Qualitatively we find the same trends also for 
ideality, Eq. || and sphericity, Eq. [7|. 

As a first possible relation between structural measures and vibration one can take the 
atomic tetragonality weighted by the amplitudes on the atoms. This would show whether 
e.g. atoms with low values of T atom participate particularly strongly in vibrations in some 
frequency range. We define 

r (") = 4E^to m i C e>)e>)> (11) 

where e n {y) stands for the three components on atom n of a vibrational eigenvector to 
frequency v and (...) denotes averaging over configurations and eigenvectors to similar 
frequencies. Taking the average value, the dashed line in Fig. [5], one observes only a slight 
variation with frequency. Only at the smallest frequencies a small upturn is found. The 
contour plot shows that the T(z/)-values fall in general into a narrow band. This is different 
for the high frequency localized modes (y > 2(e/mer) 1 / 2 ). These modes have a large spread 
of T(z/)-values without a clear preference to large or low atomic tetragonalities. This large 
spread is a direct consequence of the strong localization to even single atoms. The low 



frequency modes involve always larger numbers of atoms, from 10 upwards, and therefore 
average over many different atomic tetragonalities. 

In the previous section we noted a connection between low frequency vibrations and 
changes of structural elements. In order to quantify this notion we will now treat the 
average tetragonality, Eq. §, as a structural potential and in analogy to the usual dynamic 
matrix define a tetragonality matrix 
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Diagonalisation of this matrix gives the eigenmodes of tetragonality change and the corre- 
sponding eigenvalues, which we will denote by Gt and At, respectively. To keep in line with 
the vibrations we use a tetragonality frequency v? = vX/2%. 

In analogy to Eq. [TI|, where we defined an amplitude weighted tetragonality as function 
of frequency, we calculate the expectation value of the tetragonality matrix with respect to 
the vibrations, i.e. an amplitude weighted structural curvature, 

(ran \ 

£C(^7 e ^))- (13) 

This expectation value shows several interesting features, Fig. || Most obviously there is a 
clear more or less linear increase with frequency. This linearity breaks down at the lowest 
frequencies v < 0.2(e/ 'ma) 1 / 2 ), i.e. in the frequency range of the boson peak, where we 
find a distinct upturn. This upturn corresponds to the one in Fig. |5| but is much more 
pronounced. It clearly indicates a structural difference of the excess modes in the boson 
peak, where a small maximum resembling a boson peak is seen. It should be remembered 
that the translational invariance requires (e(z/)Te(V)) — > for v — > 0. For pure translation 
we get of course zero. Due to the limited system size sound waves below v < 0.2(e/mcr) 1 / 2 
were eliminated and we do not see the increase towards this "structural boson peak" on the 
low frequency side. The small dips of the curve for v th 0.62, 0.88, ....(e/mer) 1 / 2 coincide with 
the frequencies of the longitudinal sound waves in the SSG.Ej 

To get some deeper insight into the interplay of vibration and structure change we cal- 
culate the correlation matrix between the vibrational eigenmodes and their tetragonality 
counterparts 

(e(y)erM) = (E K^WlMf) ■ (14) 

\ na I 

The resulting correlation, Fig. |7], shows several interesting features. First there is a clear over- 
all correlation as expected from Fig. E. The correlation is highest for the highest frequency 
modes. From the participation ratiosiil on can see that both vibrational and tetragonality 
modes are localized for the highest frequencies. For the great majority of modes two groups 
can be distinguished. The largest contribution stems from a broad band stretching from the 
lowest //-values to the peak at the maximal //-values. In front of this band (higher //-values) 
there is a smaller one which can be identified as being due to longitudinal phononsE^ which 
are well separated from the other vibrations in the SSG. A third group is seen as narrow 
ridge at low v covering a major part of the v? range. This last feature shows again the 
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difference between the quasi-localized low frequency modes and the rest of the spectrum. In 
a larger system interaction will of course mix these features. This does, however, not change 
the underlying nature of the "naked" modes.o Fig. |7| does not only show separate peaks for 
the longitudinal phonons permitted by the system size but also regarding the z/^-direction 
separate "phonons" are seen, both transversal and longitudinal. Checking the participation 
ratios of ex one finds all modes with low vt to be extended, no low frequency localized modes 
are seen. The observed correlation is insufficient to predict localization at low frequencies. 
This is not too surprising as it has been observed earlierEj that these modes are produced 
by a subtle interplay of local compression and in addition a resulting soft direction in con- 
figurational space involving several atoms. The tetragonality reproduces the first feature, 
seen in the high frequency modes, but not the second one. 

To illustrate the fine details governing localization on one hand, and the stability of the 
overall correlation on the other one we repeat the above calculation for a mixed measure of 
ideality and tetragonality 
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The weighting of T and S was chosen somewhat arbitrarily to move the lowest eigenvalues to 
(2711'st) 2 ~ 0. Qualitatively the correlation, Fig. || is the same as the one for tetragonality, 
Fig. [?]. The phonons in ST-space are no longer so clearly discernible but low frequency 
localized ST-modes are found which are correlated to the low frequency vibrations. The 
occurrence of the low frequency modes for the mixed ST matrix is in agreement with 
Fig. [I] where the near zero curvature of Satom is due to changes of both tetrahedricity and 
octahedricity. The difference between Figs. [7] and |8| illustrates that the formation of low 
frequency quasi-localized modes depends much more subtly on structural details than is the 
case for high frequency modes. 

VI. CONCLUSION 

We have shown that the Voronoi-Delaunay geometrical approach gives an insight into 
the geometrical effects underlying the vibrations in the glass. The lowest frequency quasi- 
localized vibrations can be envisaged as being caused by an instability of the local geometry 
which is stabilized by the embedding lattice. A group of atoms is trapped between two 
configurations which can be considered as more perfect. We introduce different measures to 
quantify this perfectness. For the great majority of modes there is only a weak correlation 
between the amplitude on the single atoms and their perfectness. This reflects the delocali- 
sation of the modes. The high frequency localized modes which are concentrated on one or 
two atoms show a large scatter of their geometrical parameters which indicates that they 
are caused by different local distortions. At the low frequency side there is a small increase 
of tetragonality which is, however, masked by the width of the distribution. 

Introducing structural dynamic matrices correlation effects are clearly observable. These 
correlations divide the vibrations into different groups. First there are two bands of extended 
modes, longitudinal and transverse ones. The separation of these two band is due to the 
large difference in longitudinal and transverse sound velocity for the considered model. At 
high frequencies localized vibrations are correlated to high frequency structural modes. At 
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the lowest frequencies, the region of the boson peak, the vibrations show a distinctly different 
correlation behavior. This is a clear indication of their structural origin. Using a suitable 
mixture of structural measures low frequency structural modes can be defined which are 
correlated to the low frequency quasi-localized (resonant) modes. 
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FIGURES 




FIG. 1. Potential energy of atoms as function of the normal coordinate, x, corre- 

sponding to a soft mode. Solid line: average potential energy of all atoms in the system. 
(Ep t(x = 0)) ~ 7.388e. Dashed line: average potential energy of the 13 atoms most active in 
the soft mode. (E pot (x = 0))i3 « 7.364e. 
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FIG. 2. Variation of local structure in the soft mode. Black spheres: the five atoms most active 
in the mode. Gray spheres: geometrical neighbors of the active atoms. DS of perfect tetragonal 
shape (T < 0.003) are shown, a) 7 DS, b) 2 DS, c) 4 DS. 
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FIG. 3. Similar to Fig 
DS, c) 12 DS. 



DS of perfect quart-octahedral shape (O < 0.008). a) 5 DS, b) 2 
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FIG. 4. Average ideality, S^tom; of the nearest neighborhood as the function of the normal coor- 
dinate, x, of the soft mode used in Fig. (|l|). Solid line: S'atom averaged over all atoms of the system. 
(Sa,tom(x = 0)) « 0.67. Dashed line: S'atom averaged over 13 active atoms. (S a tom(£ = 0))i3 ~ 0.88. 
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FIG. 5. Amplitude weighted atomic tetragonality as function of vibrational frequency, mean 
value (dashed line) and equidistant contour lines. 
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FIG. 6. Expectation value of tetragonality with respect to vibration modes versus frequency. 
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FIG. 7. Correlation of vibrational eigenvectors, e, and the eigenvectors, ex, of the tetragonality 
matrix as function of vibrational frequency v and tetragonality frequency v^. 
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FIG. 8. Correlation as in Fig. for a mixed measure of tetragonality and ideality, Eq. |TJ 
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TABLES 

TABLE I. Values of ideality, tetragonality and sphericity for an atom in an ideal fcc-lattice, at 
the center of an icosahedron and for the glass, averaged according to Eq. ^. 





X'-'atom/ 


X-'atom/ 


(Sph atom ) 


fee lattice 


0.0 


0.039 


0.346 


icosahedron 


0.095 


0.0015 


0.325 


SS-glass 


0.66 ± 0.17 


0.023 ± 0.005 


0.358 ± 0.025 
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